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FOREWORD 


This  training  document  was  written  to  assist  users  of  Computer  Program 
HBC-2,  water  Surface  Profiles,  in  the  analysis  of  split  flows.  Financial 
assistance  for  writing  this  document  was  provided  by  the  Flood  Plain  Management 
Branch,  Office  of  the  Chief  of  Engineers. 

The  author  wishes  to  acknowledge  contributions  by  Vernon  R.  Bonner  and 
Bill  S.  Richer t  to  the  material  contained  in  this  document  and  for  their 
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APPLICATION  OP  THE  HEC-2  SPLIT  PLOW  OPTION 


INTRODUCTION 


Purpose  and  Scope. 

The  purpose  of  this  training  document  is  to  provide  information  on  the 
capabilities  and  limitations  of  the  computer  program  HBC-2  Split  Plow  option. 
This  document  provides  a  detailed  description  of  the  input  requirements,  output 
results,  computational  methods  and  example  applications  for  the  split  flow 
option.  It  is  assumed  that  the  user  has  a  knowledge  of  the  basic  HBC-2  input 

requirements.  Information  on  HBC-2  input  requirements  is  available  in  the 

\ 

"HBC-2  Hater  Surface  Profiles-Users  Manual"  (reference  1). 

\ 

\ 

Program  Documentation. 

The  split  flow  option  was  added  to  the  HBC-2  Computer  Program  by 
Modification  55  which  was  implemented  April  1982.  The  primary  documentation  for 
HBC-2  is  the  January  1981  Users  Manual  which  describes  the  program  capabilities, 
input  requirements,  and  output.  To  use  the  program,  one  would  need  a  users 
manual,  as  this  document  does  not  give  details  on  many  program  features  and 
input  formats.  The  users  manual  is  available  from  the  Hydrologic  Engineering 
Center . 
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HBC-2  SPLIT  PLOW  OPTION 


Split  flows  are  flows  that  leave  the  main  river  flow  and  take  completely 
separate  paths  from  the  path  taken  by  the  main  river  flow.  The  split  flows  may 
return  further  downstream  or  may  be  completely  lost.  The  HBC-2  program  assumes 
that  none  of  the  split  flows  return  and  automatically  reduces  the  discharges 
downstream  from  the  split  flow  location  unless  the  user  specifies  where  to 
return  the  flow  and  what  percent  of  it  to  return. 

Capabilities. 

The  HBC-2  split  flow  option  has  the  following  capabilities: 

.  Can  solve  up  to  100  separate  split  flows  simultaneously. 

.  Up  to  15  profiles  may  be  solved  in  one  execution. 

.  Three  methods  for  determining  the  split  flows  are  available. 

1.  Heir  flow  assumption 

2.  Normal  depth  assumption 

3.  Rating  curve  assumption 

.  All  or  a  percentage  of  the  split  flow  can  be  returned. 

.  Option  of  either  using  the  water  surface  or  energy 
elevation  for  computing  the  split  flows. 

.  Allows  the  use  of  rating  curves  for  starting  the  backwater. 

This  capability  is  now  a  general  HBC-2  capability  and  may 
he  used  without  having  to  use  the  split  flow  option. 


The  snllt  flow  ootion  is  compatable  with  multiple  profiles  an*!  most  of  the 


standard  HBC-2  options.  The  only  options  not  available  to  the  user,  if  the 
split  flow  option  is  used,  are  encroachment  methods  3,  4,  5,  and  6. 


There  are  many  tyoes  of  split  flows  that  occur  in  rivers.  The  following  is 
a  list  of  the  more  commonly  encountered  types  of  split  flows. 


.  Split  flows  caused  hv  islands  or  hiqh  ground. 

.  Split  flows  caused  by  levee  overtopping. 

.  Split  flows  caused  by  watershed  divide  overtopping. 
.  Split  flows  caused  by  diversion  structures. 


The  HEC-2  split  flow  option  is  capable  of  analyzing  all  of  the  above  split 
flows,  with  the  exception  of  the  island  or  high  ground  type  of  split  flows. 
Appendix  III  has  been  provided  in  this  document  to  assist  the  user  that  has  a 
split  flow  caused  by  islands  or  high  ground.  The  solution  procedure  is  the 
classic  divided  flow  analysis  technique. 
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Computation  Procedure 


The  computation  procedure  used  to  solve  the  problem  of  split  flows  is 
basically  a  one  dimensional  steady  state  method  of  trial  and  error,  as  follows: 


1.  The  program  computes  the  water  surface  Profile  and  adjusts  the 
discharges  to  reflect  the  assumed  overflows.  No  overflows  are  assumed 
to  occur  for  the  first  iteration  unless  the  user  has  specified  overflow 


2.  The  location  of  each  split  flow  reach  (Figure  1)  is  defined  in  relation 
to  the  model  by  the  upstream  and  downstream  cross  section  numbers. 
Based  on  the  cross  secion  numbers,  the  water  surface  or  energy 
elevations  are  determined  from  the  water  surface  profile  Performed  in 
step  one  and  are  used  to  calculate  the  overflow. 

3.  The  computed  and  assumed  overflows  for  each  split  flow  reach  and  their 
cumulative  values  are  compared.  If  the  difference  between  any  of  the 
values  for  computed  and  assumed  is  greater  than  2  percent,  the  program 
makes  a  new  assumption  of  overflows  and  repeats  stens  1  through  3  until 
an  acceptable  tolerance  is  met  or  unti 1  the  program  has  performed  20 


iterations 


Program  Limitations 


The  following  assumptions  are  implicit  in  the  analytical  expressions  used 
in  the  HEC-2  program  and  the  split  flow  option: 

1.  Flow  is  steady  because  time  dependent  terms  are  not  Included. 

2.  Flow  is  gradually  varied  because  a  hydrostatic  pressure  distribution  is 
assumed  in  the  energy  equation. 

3.  Flow  is  one-dimens ional  because  the  total  energy  head  is  the  same  for 
all  points  at  a  cross  section  location. 

4.  The  river  has  a  channel  slope  of  no  more  than  10  percent. 

5.  The  solit  flows  can  be  estimated  hy  the  standard  weir  equation,  the 
normal  depth  equation,  or  a  rating  curve  of  outflow  -vs-  elevation. 

6.  Submergence  of  the  overflow  weir  by  tailwater  is  in significant. 

7.  The  weir  flow  is  linearly  integrated  along  the  length  of  the  weir  based 
on  the  upstream  and  downstream  water  surface  or  energv  elevations. 

8.  The  normal  depth  conveyance  is  linearly  integrated  along  the  length  of 
the  normal  depth  cross  section  based  on  the  upstream  and  downstream 
water  surface  el >vations. 

9.  The  direction  of  the  main  stream  flow  is  at  right  angle  to  the  split 
flow. 

10.  Split  flow  is  controlled  by  either  the  water  surface  elevation  or  the 
energy  elevation. 

11.  Flow  boundaries  are  fixed.  This  is  to  say  that  the  cross  section  and 
weir  geometries  do  not  erode  or  change  with  time. 
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Split  flows  that  are.  not  returned  to  the  system  are  removed  from  the 
entire  model  downstream  from  the  split  flow  location.  because  HEC-2 
does  not  take  into  account  the  variation  of  flow  with  time,  the  user 
must  be  careful  in  using  the  HEC-2  split  flow  ootion  in  cases  where  the 
split  flow .  at  an  upstream  location  does  not  have  a  constant  effect  on 
the  peak  discharges  further  downstream.  An  explanation  of  this  problem 
and  a  solution  strategy  is  given  in  the  section  titled  Hvdrologic 


Considerations 


HEC-2  SPMT  PLOT  INPOT  DESCRIPTION 


The  solit  flow  innut  data  must  be  enters'1  as  the  first  data  of  an  HBC-2 
run,  with  the  exception  of  the  ED  card  (B0IT2  nroqram) .  The  split  flow  innut 
data  must  alwavs  start  with  an  SF  card  and  end  v.  th  an  EE  card.  The  EE  card 
must  be  followed  by  the  reoular  HEC-2  innut  data  cards. 

The  split  flow  input  data  uses  the  standa  d  hBC  innut  format,  in  which  the 

first  field  contains  a  two  character  car  •  v  i  ^  i  *=»r #  and  has  six  columns  for 

•  data  while  the  next  nine  fields  contair  -•  -  each.  The  above  format  is 

for  numeric  data.  The  format  for  title  cnr  •-  •  • -->  columns  for  card  identifier 

and  the  remaining  seventy  eiqht  columns  for  aloha-numeric  information. 

The  split  flow  data  input  varies  from  the  standard  HEC-2  format  in  that  a 
set  of  data  cards  is  always  preceded  bv  a  title  card  and  the  title  card  is 
required.  The  order  in  which  the  cards  are  entered  must  be  followed  exactly  for 
each  set  or  qroup  of  cards.  The  sets  of  cards  need  not  be  in  any  specific 
order,  but  it  is  recommended  that  split  flow  reaches  be  entered  in  a  downstream 
direction  to  make  the  output  more  readable.  A  detailed  description  of  the  split 
flow  data  card  inout  format  is  qiven  in  Appendix  I.  The  following  is  a 
description  of  the  qeneral  card  types. 

*■  Split  Plow  Title  Card. 

The  SF  card  is  a  title  card  used  to  activate  the  split  flow  option.  It  is 

a  required  card  and  must  be  the  first  input  data  card  in  the  deck. 
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Job  Card  Set 


The  job  card  set  consists  of  the  JC  and  jp  cards.  The  job  card  set  is 
optional  and  is  used  to  control  the  processing  of  the  split  flow  data.  A  job 
card  set  may  be  pladed  anywhere  in  the  sollt  flow  Input  data.  The  JC  card  is 
the  first  card  of  the  set  and  is  used  as  a  title  card.  The  JP  card  is  used  to 
set  the  level  of  printout,  allowed  error  tolerance,  maximum  number  of 
iterations,  use  of  either  water  surface  or  energy  elevation,  and  percent  of 
split  flow  to  return. 

Weir  Reach  Card  Set. 

The  weir  reach  set  is  composed  of  three  types  of  cards.  The  first  card  of 
the  set  is  the  TW  card.  The  TW  card  is  a  title  card  and  must  be  the  first  card 
of  the  set.  The  second  card  is  the  WS  card,  which  contains  information  dealing 
with  the  number  of  points  describing  the  weir,  weirflow  coefficient,  location  of 
the  downstream  and  upstream  limits  of  the  weir  in  relation  to  section  numbers 
and  the  section  number  where  the  flow  returns.  The  third  card  of  the  set,  the 
WC  card,  describes  the  weir  geometry  by  the  use  of  station  and  elevation 
coordinates.  The  coordinate  points  must  start  at  the  downstream  end  and  proceed 
in  an  upstream  direction. 
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card  of  the  set  is  the  TN  card,  which  is  a  title  card  and  must  be  the  first  card 
of  the  set.  The  second  card,  the  NS  card,  is  similar  to  the  WS  card,  with  the 
exception  that  instead  of  havinq  the  wetrflow  coefficient,  it  has  the  enerqy 
slope  and  'n'  value.  The  third  card,  the  NG  card,  is  used  to  describe  the 
normal  depth  cross  section  aeometry  by  the  use  of  station  and  elevation 
coordinates,  with  the  coordinates  starting  at  the  downstream  end. 

Diversion  Reach  Card  Set. 

The  diversion  reach  set  is  composed  of  three  types  of  cards.  The  first 
card  of  the  set  is  the  TC  card,  which  is  a  title  card  and  mist  be  the  first  card 
of  the  set.  The  second  card,  the  CS  card,  contains  information  dealing  with  the 
number  of  points  describing  the  diversion  rating  curve,  loactlon  of  the 
diversion  in  relation  to  section  numbers,  and  the  section  number  wN»re  the 
diverted  flow  returns.  The  third  card,  the  CR  card,  is  used  to  def,crib»  the 
ratinq  curve  bv  the  use  of  discharge  and  elevation  coordinates. 

End  Split  Flow. 

The  EE  card  is  required  to  terminate  the  split  flow  input  data. 


Additional  HEC-2  Cards 


Several  additional  input  data  cards  have  been  added  to  the  standard  HBC-2 
input  to  facilitate  the  use  of  the  split  flow  ootion.  The  new  cards  are  the  JR, 
JS,  and  RC  cards. 

The  JR  card  is  used  to  input  a  ratine  curve  that  is  used  to  start  the 
backwater .  The  JR  card  follows  the  Jl  card  and  is  read  when  the  Jl  card 
variable  STRT  is  greater  than  one.  The  STRT  value  in  this  case  is  used  to 
indicate  the  number  of  rating  curve  values  that  will  be  read  on  the  JR  card. 

The  RC  card  is  used  to  input  a  rating  curve  at  any  cross  section,  which 
will  be  used  instead  of  a  calculated  backwater  answer.  It  operates  in  the  same 
manner  as  an  X5  card. 

The  JS  card  is  used  to  specify  the  starting  assumed  split  flow  for  each 
reach  defined  in  the  split  flow  data  sets.  The  JS  card  follows  the  Jl  or  JR 
cards.  It  is  an  optional  card  and  if  omitted,  the  program  assumes  for  its  first 
trial  that  no  flow  is  being  lost. 


GENERAL  MODELING  CONSIDERATIONS 


The  snlit  flow  option  requires  that  the  cross  section  numhers  (Xl  card) 
continually  increase  in  value  from  downstream  to  uostream.  It  is  recommended 
that  station  values  in  feet  or  miles  measured  along  the  main  channel  he  used  for 
cross  section  numbers.  This  is  imoortant  because  the  HBC-2  proqram  uses  the 
cross  section  numbers  specified  on  the  split  flow  reaches  to  determine  what 
water  surface  or  energy  qrade  line  elevations  to  use  in  determining  the  amount 
of  flow  lost.  Each  split  flow  reach  is  located  based  on  cross  section  numhers 
at  the  downstream  and  uostream  ends  of  the  reach.  If  the  downstream  and 
upstream  split  flow  locations  do  not  match  a  cross  section  number,  the  proqram 
will  linearly  interpolate  between  cross  section  numbers  to  determine  the  water 
surfaces  or  energy  grade  line  elevations.  A  wise  practice  would  be  to  start  and 
end  the  split  flow  reaches  at  cross  section  numbers  that  appear  in  the  NBC-2 
hydraulic  model. 

Split  Plow  Reach  Lenqth  Consideration. 

The  overflow  reaches  should  be  kept  as  short  as  possible.  The  longer  the 
split  flow  reaches,  the  less  accurate  will  be  the  split  flows  and  backwater 
answers.  The  split  flow  problem  can  be  compared  to  the  integration  of  a  curve 
by  the  Trapesoidal  Rule,  which  is  an  approximate  method  of  subdividing  the  curve 
into  a  number  of  straight  line  segments  and  calculating  the  areas  under  each 
straiqht  line  segment  and  then  adding  up  the  subareas  to  calculate  the  total 
area  under  the  curve.  The  smaller  the  straight  line  segments  and  the  larger  the 
number  of  segments,  the  more  accurate  will  be  the  calculated  area.  The  same 
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RIVER  DISTANCE 


c  loqic  a op lies  to  the  calculation  of  the  solit  flows.  The  shorter  the  split  flow 

! 

reaches  and  the  additional  cross  sections  defining  them,  the  more  accurate  will 
be  the  calculated  outflows  and  water  surface  elevations. 

The  effect  of  subdividing  a  split  flow  reach  and  addition  of  cross  sections 
can  be  seen  in  the  following  examole.  Figure  2  shows  a  profile  Plot  of  a  stream 
and  a  split  flow  reach  which  has  a  total  reach  length  of  1000  feet.  The  split 
flow  reach  was  modeled  in  four  different  ways.  The  first  model  represented  the 
split  flow  reach  as  a  single  1000  foot  long  reach.  The  second  model  divided  the 
split  flow  reach  into  two  500  foot  reaches.  The  third  r.odel  divided  the  split 
flow  reach  into  four  250  foot  reaches  and  the  fourth  model  divided  it  into  ten 
100  foot  reaches.  Note  that  each  reach  has  a  crons  section  defininq  it  upstream 
and  downstream.  Additional  cross  sections  must  he  used  because  the  split  flow 
discharges  are  calculated  based  on  a  linear  interpolation  of  the  upstream  and 
downstream  cross  section  water  surface  or  energy  elevations.  The  effect  of 
dividing  up  the  split  flow  reach  into  shorter  and  shorter  segments,  is  to 
produce  answers  which  are  more  accurate.  The  results  from  these  four  models  are 
presented  on  Fiaure  3,  which  shows  plots  of  the  flow  in  the  main  river  -vs- 
river  distance.  The  results  clearly  show  that  the  1000  foot  reach  should  be 
divided  into  shorter  reaches.  The  fourth  model  produced  the  most  accurate 
results,  but  it  is  also  clear  from  the  plotted  results,  that  the  second  and 
third  models  gave  acceptable  results. 

The  example  also  shows  that  the  upstream  portion  of  a  uniform  split  flow 
reach  should  he  divided  into  shorter  segments  than  the  downstream  portion, 
because  a  larger  proportion  of  the  flow  will  be  lost  on  the  upper  portion. 
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Hydroloqic  Considerations. 

A  common  practice  for  hydraulic  analysis  is  to  use  the  flood  hydroqraph 
peak  discharge  values  for  the  backwater  calculations  and  ignore  the  fact  that  in 
most  cases  the  peaks  of  the  hydroqraphs  do  not  occur  at  the  same  time.  In  many 
cases  the  effect  of  this  assumption  on  the  calculate*  water  surface  profile  is 
well  within  the  accuracy  for  profile  calculations.  Most  HT5C-2  models  are 
assembled  in  this  manner.  When  it  comes  to  a  solit  flow  analysis,  the  effect  of 
using  the  peak  djgcharaes  as  steady  flow  should  not  he  iqnored. 

The  HBC-2  split  flow  option,  as  a  default,  reduces  all  discharge  values 
downstream  from  a  split  flow  reach  by  the  calculated  split  flow  values.  For 
example,  if  1000  cfs  is  lost  at  the  headwaters  of  a  stream,  the  program  will 
reduce  all  the  discharges  by  1000  cfs  all  the  wav  downstream  to  the  start  of 
model  or  to  the  section  number  where  the  user  has  specified  that  it  returns. 
The  onlv  way  to  determine  the  validity  of  reducing  the  discharges  in  this  manner 
is  to  look  at  the  entire  hydrograph,  to  see  how  the  loss  from  the  split  flow 
effects  the  Peak  discharges  downstream.  It  may  be  that  the  peak  discharge 
further  downstream  will  onlv  be  partially  effected  by  the  split  flow  loss.  The 
peak  discharges  downstream  may  be  more  dependent  on  the  timing  of  tributary  and 
local  inflows. 

To  illustrate  this  problem,  the  following  simple  example  is  presented.  As 
shown  on  Figure  4,  the  Upper  Main  Stem  has  a  split  fow  reach  just  upstream  from 
the  confluence  with  Tributary  No  1.  Fiqure  5  is  a  plot  of  the  original 
hydrographs  for  the  Upper  Main  Stem,  Tributary  No  l,  and  their  combined 
hydrographs  on  the  Lower  Main  Stem.  If  flows  that  exceed  3000  cfs  on  the  Upper 
Main  Stem  are  lost,  then  the  modified  hydrographs  plotted  on  Figure  6  occur. 
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Cowparing  the  original  hydrographs  with  the  modified  shows  the  peak  discharge  on 


the  Lower  Main  Stem  reduced  by  1000  cfs,  peak  flow  on  Tributary  No  1  unaffected, 
and  the  Peak  discharge  on  the  Lower  Main  Stem  reduced  by  S00  cfs.  The  HEC-2 
split  flow  option  woi/ld  have  reduced  the  Lower  Main  Stem  peak  discharge  by  1000 
cfs  and  not  by  the  correct  mount  of  500  cfs. 


FIGURE  4 

Consideration  must  also  he  qiven  to  the  effect  of  storage  routinqs.  The 
effect  of  routing  a  flood  wave  Is  to  reduce  the  peak  discharge.  Therefore  a 
peak  discharge  further  downstream  based  on  routinq  will  not  he  reduced  bv  the 
total  amount  that  was  lost  upstream. 
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The  following  procedure  can  be  used  to  account  for  the  hydroloqic  aspects 
of  a  split  flow  problem. 

1.  Hake  an  initial  HBC-2  run  which  reduces  the  split  flows  all  the  wav 
downstream. 

2.  Alter  the  hydroloqic  model  (REC-1) ,  at  the  split  flow  reach,  to  reflect 
the  lost  flow  and  execute  it  again.  Analyze  the  effect  that  the  split 
flows  had  on  the  peak  discharges  downstream  from  the  split  flow  reach. 

3.  If  the  assumption  of  reducinq  the  lost  flow  is  valid,  then  no  further 
analysis  are  needed.  If  the  lost  flow  is  only  partially  effective 
downstream,  then  the  lost  flows  should  be  returned  further  downstream. 
This  can  be  accomplished  by  modifying  the  HBC-2  discharge  cards  to 
reflect  the  expected  reductions  downstream  from  the  split  flow  reach  and 
returning  the  split  flows  back  into  the  model. 

4.  An  HBC-2  run  should  be  made  and  steps  2  and  3  repeated  until  an 
acceptable  solution  is  achieved. 

This  procedure  is  applicable  for  simple  split  flow  problems  which 
do  not  have  more  than  three  or  four  separate  split  flow  reaches.  For 
the  more  complicated  split  flow  problem,  an  unsteady  state  oroqram 
should  be  used.  The  DWOPBR  (Dynamic  Have  Operational  Model)  program 
developed  by  the  National  Heather  Service  (reference  8)  is  a  one 
dimension  unsteady  state  program  that  can  be  used  to  solve  the  split 
flow  problem. 
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APPENDIX  I 


SPLIT  FLOW  INPUT  AND  OUTPUT  DATA  DESCRIPTION 


ft *?«¥!»  V-' 


SLIT  FLOW  TITLE  CARD 


CARD  SF  -  SPLIT  FLOW  CARD  (REOtJTRED  IF  SPLIT  FLOW  OPTION  IS  TO  BE  DSHD) 

The  SF  card  is  used  to  flap  the  split  flow  option.  Only  one  SF  card  can 
be  used.  The  SF  card  has  to  be  the  first  card  in  an  HEC-2  deck. 


Field  Variable  Value 


Descr iption 


0 


IA 


SF  Card  identification  characters 


1-10 


Aloha-numeric  Title  data. 


CARD  JC  -  TITLE  JOB  CARD  FOR  SPLIT  FLOW  (OPTIONAL) 

The  JC  card  is  used  to  indicate  that  JP  card  follows.  The  JP  card  must 
follow  the  JC  card. 

Field  Variable  Value  Description 

0  IA  JC  Card  identification  characters 

1-10  Alpha-numeric  Title  data. 


t 

* 

i 
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CARD  JP  -  JOB  PARAMETER  CARD 


The  JP  card  is  used  to  set  several  job  parameters  dealing  with  the 
split  flow  computations.  The  JC  and  JP  cards  are  optional  and  can  he 
placed  anywhere  in  the  split  flow  data  or  completely  left  out.  Thev 
should  be  placed  normally  after  the  SF  cards. 


Field 

Variable 

Value 

Description 

0 

IA 

JP 

Card  identification  characters. 

1 

1SFTR 

0 

Printout  control  of  split  flow 
computations  will  be  held  to  a  minimum. 

1 

Trace  each  snl.it  flow  iteration. 

10 

Trace  both  the  split  flow  and  backwater 
iterations. 

2 

AEROR 

0 

The  proqram  will  use  a  value  of 

2  percent  allowed  error  for  convergence. 

+ 

The  user  may  specify  the  allowed  percent 
tolerance  for  convergence. 

3 

NAITER 

0 

The  maximum  number  of  iterations 
for  split  flow  to  be  executed  per 
profile  (20  is  the  default  value). 

+ 

The  user  may  specify  the  maximum 
number  of  iterations. 

4 

IUEG 

-1,0 

The  proqram  will  use  the  water 
surface  to  determine  the  overflow. 

1 

The  oroqram  will  use  the  energy  grade 
line  to  determine  the  overflow. 

5 

PERFR 

0 

One  hundred  percent  of  the  overflow 

is  to  be  returned  at  SNOFR  (WS.4,  NS. 4, 
and  CS.4) . 


+  Percent  of  overflow  to  be  returned  at 

SNOFR  (WS.4,  NS. 4,  and  CS.4). 
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CARD  TW  -  TITLF  CARD  FOR  WEIR  DOCATION 


The  TW  card  is  required  for  each  set  of  weir  outflow  data  set.  The 
TW  card  must  he  followed  hy  a  set  of  WS  and  WC  cards. 


Field 

Variable 

Value 

Description 

0 

IA 

TW 

Card  identification  characters 

1-10 

Alpha-numeric  Title  data. 
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CARO  WS  -  WEIR  PARAMETER  DATA  CARD 


The  WS  car d  is  required  for  each  TW  card  used  and  must  follow  it. 

The  WS  card  contains  information  dealinq  with  the  number  of  points 
describinq  the  weir,  weir  flow  coefficient,  location  of  the  upstream 
and  downstream  limits  of  the  weir  in  relation  to  section  numbers  as 
used  in  the  Xl  cards,  and  the  section  number  where  the  flow  returns. 
If  the  flow  does  not  return,  a  value  of  -1  should  he  used.  It  is 
required  that  the  section  numbers  used  to  set-up  the  backwater 
model  increase  from  downstream  to  uostream.  The  same  rule  applies 
for  supercritical  models. 


Pield 

Variable 

Value 

Description 

0 

IA 

WS 

Card  identification  characters. 

1 

NWPL 

+ 

Number  of  coordinate  points  that  describe 
the  weir  on  the  WC  card. 

2 

nSSNO 

o,+ 

Downstream  section  number  where  the 
first  weir  coordinate  apolies. 

3 

USSNO 

0,+ 

Upstream  section  number  where  the 
last  weir  coordinate  applies. 

4 

SNOFR 

0,+ 

Section  number  where  the  lost,  weir 
flow  returns. 

-1 

The  weir  flow  does  not  return. 

5 

COEFL 

+ 

Coefficient  of  discharge  Cor  use  in  weir 
flow  equation. 

6-10 

Not  used. 
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CARD  WC  -  THE  WEIR  COORDINATE  CARD 


The  WC  card  is  used  to  input  the  weir  coordinates.  The  weir  coordinates 
must  start  at  the  downstream  end  and  proceed  upstream.  The  maximum  number 
of  coordinates  is  100. 


Field 

Var table 

Value 

Description 

0 

IA 

WC 

Card  identification  characters. 

1,3,5, 

7,9 

STA(I) 

+ 

Station  Value  of  weir  coordinate. 

2,4,6, 

8,10 

ELO(I) 

+ 

Elevation  value  of  weir  coordinate. 

CARD  TN  -  TITLB  CARD  FOR  NORMAL  DEPTH  LOCATION 

The  TN  card  is  required  for  each  set  of  normal  depth  outflow  data  set. 
The  TN  card  must  be  followed  by  a  set  of  NS  and  NO  cards. 


Field 

Variable 

Value 

Descr iDtion 

0 

IA 

TN 

Card  identification  characters 

1-10 

Alpha-numeric  Title  data. 

<  .tva-  «*-!* 


CARP  NS  -  NORMAD  DEPTH  PARAMETER  DATA  CARD 


The  MS  card  is  similar  to  the  WS  card  with  the  exceotlon  that  instead  of 
having  the  weir  flow  coefficient,  it  has  the  energy  slope  and  ‘n*  value. 


Pi  eld 

Variable 

Value 

0 

IA 

NS 

1 

NWPL 

+ 

2 

OSS  NO 

o,+ 

3 

USSNO 

0,+ 

4 

SNOPR 

0,+ 

-1 

5 

XNVND 

+ 

6 

SLOPND 

+ 

7-10 


Description 

Card  identification  characters. 

Number  of  coordinate  points  that  describe  the 
normal  depth  flow  cross  section  on  the  NG  card 

Downstream  section  number  where  the  first 
coordinate  point  on  the  NG  card  applies. 

Upstream  section  number  where  the  last 
coordinate  point  on  the  NG  card  applies. 

Section  number  where  the  lost  flow  returns. 

The  lost  flow  does  not  return. 

The  'n'  value  to  be  used  for  normal  depth 
calculation. 

The  energy  slope  to  he  used  for  normal 
deoth  calculations. 

Not  used. 
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CARD  NO  -  THE  GROUND  COORDINATE  CARD 


The  NG  card  is  used  to  input  the  normal  depth  cross  section  coordinates. 
The  coordinate  must  start  at  the  downstream  end  and  proceed  upstream. 

The  maximum  number  of  coordinates  is  100. 


field 

Var iable 

Value 

0 

IA 

NG 

1#3,5, 

7,9 

STA(I) 

+ 

2,4,«, 

8,10 

ELO(I) 

♦ 

Descr iptlon 

Card  identification  characters. 
Station  Value  of  cross  section. 

Elevation  value  of  cross  section 
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CARP  TC  -  TITLE  CARO  FOR  RATING  CURVE  U3CATION 

The  TC  card  is  required  for  each  set  of  ratinq  curve  outflow  data 
set.  The  TO  card  must  be  followed  by  a  set  of  CS  and  CR  cards. 


Field 

Variable 

Value 

Descr iption 

0 

IA 

TC 

Card  identification  characters 

1-10 

Aloha-numeric  Title  data. 

CARD  CS  -  RATING  CURVE  PARAMETER  DATA  CARD 

The  CS  card  is  similar  to  the  WS  card  with  the  exception  that  the 
location  (unstream  and  downstream)  is  a  nolnt  location  and  therefore 
the  value  entered  for  WSSNO  and  OSSNO  should  normally  be  equal. 


Field 

Variable 

Value 

Description 

0 

IA 

CS 

Card  identification  characters. 

1 

MWPL 

+ 

Number  of  discharqe  elevation  pairs  to  be 
read  from  the  CR  cards  to  follow. 

2 

DSSNO 

o,+ 

Downstream  section  number  where  the 
ratinq  curve  applies. 

3 

USSNO 

0,+ 

Upstream  section  number  where  the 
ratinq  curve  applies. 

4 

SNOFR 

+ 

*  H 

O  1 

Section  number  where  the  lost  flow  returns. 

The  lost  flow  does  not  return. 

5-10  Not  qsed. 
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CARO  CR  -  RATING  CURVE  CARO 


The  CR  card  is  used  to  input  the  rating  curve  of  outflows.  The 
location  of  the  rating  curve  ha*  to  he  at  a  specific  location  on 
the  river.  Therefore  the  location  has  to  he  specified  at  only  one 
point.  The  variables  DSSNO  and  USSNO  should  he  set  eoual.  If  they 
are  not,  the  program  will  use  the  mean  of  the  two  locations. 

The  maximum  number  of  rating  curve  points  is  100. 


field 

Variable 

Value 

Description 

0 

IA 

CR 

Card  identification  characters. 

1,3,5, 

7,9 

STA(I) 

+ 

Discharge  values  for 

rating  curve 

2,4,6, 

8,10 

ELO(I) 

+ 

Elevation  values  for 

rating  curve 

CARD  EE  -  END 

OP  SPLIT  PLOW 

DATA  CARD 

The  EE 

data. 

card. 

card  is  required  to  terminate  the 
The  EE  card  should  he  in  front  of 
such  as  the  AC,  C,  or  T1  cards. 

reading  of  the  split  flow 
the  first  regular  HEC-2 

rieid 

Var iable 

Value 

Descriotion 

0 

IA 

EE 

Card  identification  characters. 

1-10 

Not  used. 
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MODIFICATION  TO  J1  CARO 


The  J1  card  variable  STRT  (Field  5)  has  been  altered  so  that  the  program 
will  accept  a  rating  curve  as  a  startinq  backwater  condition.  The  option 
is  activated  by  entering  in  field  five  (STRT  variable)  of  the  .71  card  the 
number  of  discharge  elevation  rating  curve  points.  The  rating  curve  is 
entered  after  the  Jl  card  using  JR  cards. 


CARO  JR  -  STARTING  RATING  CURVE  CARO 


The  JR  cards  are  used  to  input  a  starting  rating  curve.  A  set  can  be 
placed  for  each  orofile  being  run.  Thev  must  follow  the  Jl  card  and 
the  number  of  rating  curve  points  must  be  greater  than  two.  It  is 
required  that  the  number  of  rating  curve  points  be  entered  on  the  Jl 
card,  field  five.  A  maximum  of  twenty  discharge  elevation  values 
is  allowed.  The  program  linearlv  interpolates  between  given  ratinq 
curve  values  and  extrapolates  for  values  outside  the  rating  curve. 


Field 

Var  table 

Value 

Description 

0 

IA 

JR 

Card  identification  characters 

1,3,5, 

7,9 

ojl(i) 

+ 

Discharge  values. 

2,4,6, 

8,10 

XJ1(I) 

-,0,+ 

Water  surface  elevation  values 

1-11 


1 


JS  -  CARD  STARTING  SPLIT  PLOW  ASSUMPTION  CARD 

The  JS  car*  is  used  to  specify  the  startinq  assumed  lost  discharqes 
for  each  reach  defined  in  the  split  flow  data  set.  If  the  JS  card 
is  not  entered  for  a  profile,  then  the  proqraia  assumes  that  the  first 
trial  assumed  lost  flow  is  zero  for  all  the  split  flow  reaches.  The 
JS  card  should  follow  the  J1  card  or  the  JR  card  if  used.  A  maximum 
of  100  values  are  allowed. 


Field 

Variable 

Value 

Description 

0 

TA 

JS 

Card  identification  characters. 

1 

N 

+ 

Number  of  assumed  lost  discharqes 

to  read 

2 

ARL0(4,1) 

+ 

Assumed  lost  discharge  for  first  : 

reach. 

3 

• 

• 

ARL?>  (4,2) 

• 

e 

+ 

• 

e 

Assumed  lost  discharge  for  second 

e 

• 

reach . 

• 

• 

ARL0<4,N) 

• 

♦ 

• 

Assumed  lost  discharge  for  last 
reach. 

Continue  on 

in  field  one  of 

additional 

JS  cards  up  to  ARI£)(4,N). 
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CARD  RC  -  RATING  CURVE  CARD 


The  RC  card  can  be  entered  at  anv  cross  section  and  the  program  will 
determine  the  water  surface  elevation  based  on  the  rating  curve  and 
not  on  backwater  computations.  The  RC  card  should  be  placed  after  the 
Xl  card.  A  maximum  of  twenty  discharge  elevation  values  are  allowed. 
The  program  linearly  interpolates  between  given  rating  curve  values 
and  extrapolates  for  values  outside  the  rating  curve. 


Field 

Variable 

Value 

0 

IA 

RC 

1 

NRCP 

+ 

2 

PRC(1) 

+ 

3 

XRC(l) 

-#o,+ 

4 

ORC (2) 

♦ 

5 

XPC{2> 

-,0,+ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

e 

ORC (NRCP) 

XRC (NRCP) 

• 

♦ 

-,o,+ 

Description 

Card  identification  characters. 

Number  of  rating  curve  points  being 
read  in. 

Discharge  value. 

Water  surface  elevation  value. 
Discharge  value. 

Water  surface  elevation  value. 


Last  discharge  value. 

Last  water  surface  elevation  value. 


Continue  on  in  field  one  of  additional  RC  cards  up  to  ORC (NRCP)  and 
XRC(NRCP) . 


HEC-2  SPLIT  FLOW  OPTION  OUTPUT  VARIABLES 


Var  table 


Descr lot ton 


ASQ 

O  COMP 

ERRAC 

TASQ 

TOO 

TABER 

NITER 

OSWS 
tlSWS 
DSSNO 
USSNO 
TOTAL  AREA 

AVG  VELOCITY 
MAX  DEPTH 

AV  DEPTH 

TOP  WIDTH 
TOP  WIDTH 


The  assume*  spilt  flow  value  use*  by  the  program  to  compute  the 
water  surface  elevations. 

The  compute*  split  flow  value  base*  on  the  compute*  water  surface 
elevation. 

The  percent  error  between  the  assume*  *lscharge  an*  comoute* 
discharge. 

The  total  assume*  split  flow  for  the  entire  HEC-2  model. 

The  total  compute*  split  flow  for  the  entire  HEC-2  model. 

Percent  of  error  between  the  total  assume*  split  flow  an*  total 
compute*  split  flow. 

The  number  of  iterations  that  the  program  has  execute*  in  computing 
the  answer. 

The  compute*  downstream  water  surface  elevation. 

The  compute*  upstream  water  surface  elevation. 

The  downstream  section  number  where  the  split  flow  reach  begins. 

The  upstream  section  number  where  the  split  flow  reach  ends. 

The  total  cross  sectional  area  for  a  normal  depth  overflow 
reach . 

The  averaqe  velocity  of  the  normal  depth  overflow  reach. 

The  maximum  depth  that  occurs  on  the  normal  depth  overflow 
section. 

The  average  deoth  of  flow  for  the  normal  depth  section  base* 
on  the  total  area  divide*  by  the  water  surface  topwidth. 

The  width  of  the  normal  depth  overflow  section. 

The  width  of  the  overflow  section  base*  on  the  compute* 
water  surface. 
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APPENDIX  II 

SPLIT  FLOW  EXAMPLE  PROBLEM 


EXAMPLE  OF  INPUT  PREPERATION 


The  following  problem  Is  provided  to  illustrate  the  input  preoaration 
repaired  when  using  the  solit  flow  option.  The  input  is  shown  on  Figure  9  and 
is  described  below.  A  comolete  HEC-2  computer  outout  listing  for  this  example 
is  also  provided. 

A  plan  view  of  the  levee  system  and  floodway  of  the  Red  Fox  River  is  shown 
in  Figure  7.  A  profile  view  of  the  stream  bed,  levees  and  overflow  weir  are 
shown  in  Figure  8.  As  can  be  seen  on  the  profile  view,  the  right  bank  or  south 
levee  is  the  critical  levee  for  a  solit  flow  analysis.  This  is  because  it  has 
the  overflow  weir  and  is  several  feet  lower  than  the  north  levee. 

The  starting  water  surface  for  this  example  will  be  based  on  a  normal  depth 
calculation  using  a  slope  of  .005  ft/ft.  The  weir  coefficient  for  the  levee 
will  be  3.4  and  for  the  overflow  weir,  2.7.  The  weir  coefficient  for  the 
overflow  weir  is  low  to  account  for  the  submergence  caused  by  the  tailwater  in 
the  flood wav. 

The  first  card  of  the  split  flow  inout  data  is  the  SF  card.  The  split  flow 
reach  in  this  example  is  divided  into  three  shorter  reaches.  The  first  solit 
flow  reach  to  be  modeled  is  the  most  upstream  reach  which  lies  between  sections 
3  and  4.  The  TW  card  follows  the  SF  card  and  is  used  to  identify  the  reach. 
The  WS  card  follows  the  TW  card  and  is  used  to  specify  the  number  of  coordinates 
that  describe  the  levees  geometry,  location  of  the  downstream  and  upstream 
limits  and  the  section  number  where  the  solit  flow  returns.  The  solit  flow  does 
not  return  and  a  value  of  -1  is  used  to  so  indicate.  Station  and  elevation 
coordinate  data  is  entered  on  the  WC  card  to  describe  the  weir  geometrv.  The 


coordinate  data  start  at  section  3  and  proceed  upstream  to  section  4 


The  same  procedure  used  for  the  first  upstream  reach  is  used  for  the 
overflow  weir  reach,  which  lies  between  sections  2  and  3,  and  for  the  downstream 
levee  reach  which  lies  between  section  1  and  2. 

In  this  example  problem  it  was  assumed  that  the  weir  equation  was  the  more 
anorooriate  to  use  for  calcul.atinq  the  split  flows.  If  a  normal  deoth  or 
diversion  assumption  had  been  preferred,  then  the  onlv  difference  in  input  would 
have  been  that  the  normal  depth  or  diversion  split  flow  cards  would  have  been 
used. 
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APPENDIX  III 


GRAPHICAL  METHOD  FOR  SOLVING  ISLAND  DIVIDED  FLOWS 


ISLAND  SPLIT  FLOW 


Where  an  island  or  other  obstruction  in  the  river  separates  flow  into 
two  or  more  channels  over  a  substantial  length  (several  cross  sections) , 
the  quantity  of  water  passing  on  each  side  of  the  island  must  be  determined 
since  total  energy  loss,  past  the  island,  must  be  the  same  for  both  sides. 

The  example  in  fig.  III.l  illustrates  how  to  solve  the  divided  flow 
problem  graphically.  The  total  discharge  is  proportioned  between  the  north 
and  south  channels,  arbitrarily.  The  water  surface  elevation  for  the  total 
flow  is  determined  for  river  mile  10.0,  and  a  water  surface  profile  is 
calculated  for  each  assumed  discharge  through  the  north  channel  and  through 
the  south  channel.  The  resulting  potential  water  surface  elevations  at  river 
mile  10.8  are  plotted  in  fig.  III. 2.  A  "total"  discharge  curve  is  obtained 
at  river  mile  10.8  by  summing  north  and  south  channel  discharges  for  common 
water  surface  elevations.  This  total  flow  curve  intersects  the  total  river 
discharge,  5000  cfs,  at  elevation  104.32,  thereby  defining  the  upstream 
water  surface  elevation.  By  intersecting  the  north  and  south  channel  curves 
at  elevation  104.32,  their  respective  discharges  can  be  read  from  the  figure. 


Using  the  data  for  the  discharge  of  5,000  units  tabulated  below,  determine  the  water- 
surface  elevation  at  river  mile  10.8  when  the  water-surface  elevation  at  river  mile 
10.0  is  100.0.  What  is  the  discharge  in  the  channel  north  of  the  island? 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whit  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEPORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

Tra -firing  TVirnmenf  Nn.  18 

s.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

Application  of  the  HEC-2  Split  Flow  Option 

S.  TYPE  OF  REPORT  6  PERIOO  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfoJ 

Alfredo  E.  Montalvo 

S.  CONTRACT  OR  GRANT  NUMBERS 

•  ■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Corps  of  Engineers 

The  Hydrologic  Engineering  Center 

609  Second  Street,  Davis,  CA  95616 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

13.  NUMBER  OF  PAGES 

Sfi 

14.  MONITORING  AGENCY  NAME  6  AOORESSflt  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (of  thle  report) 

Unclassified 

fsa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (o I  thit  Report) 


Distribution  of  this  publication  is  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  l/io  abtlrmct  entered  In  Block  30,  II  different  from  Report) 


16.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  tovetam  aide  i/n*c*iu/y  and  Identity  by  block  number) 

Digital  Computer,  Split  Flows,  HEC-2,  Backwater,  Open-Channel  Hydraulics, 
Lateral  Weir  Overflows 


20.  ABSTRACT  (Continue  on  reveree  tide  If  neceeemty  end  Identlfr  bp  block  mmber)  - 

Training  document  designed  to  assist  the  user  of  computer  program  HEC-2 
split  flow  option  and  to  acquaint  them  with  its  capabilities  and  limitations. 
Provide  a  detailed  description  of  the  HEC-2  split  flow  input  requirements, 
output  results,  computation  method  and  example  applications. 


DO  i'SJTt.  1473  EDITION  OF  1  NOV  6S  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  DM*  Entered) 


